Quasirelativistic quasilocal finite wave-function collapse model 



I. INTRODUCTION 



Philip PearlCl 

Department of Physics, Hamilton College, Clinton, NY 13323 
(Dated: February 1, 2008) 

A Markovian wave function collapse model is presented where the collapse-inducing operator, con- 
structed from quantum fields, is a manifestly covariant generalization of the mass density operator 
utilized in the nonrelativistic Continuous Spontaneous Localization (CSL) wave function collapse 
model. However, the model is not Lorentz invariant because two such operators do not commute at 
spacelike separation, i.e., the time-ordering operation in one Lorentz frame, the "preferred" frame, 
is not the time-ordering operation in another frame. However, the characteristic spacelike distance 
over which the commutator decays is the particle's Compton wavelength so, since the commutator 
rapidly gets quite small, the model is "almost" relativistic. This "QRCSL" model is completely 
, finite: unlike previous, relativistic, models, it has no (infinite) energy production from the vacuum 

state. QRCSL calculations are given of the collapse rate for a single free particle in a superposition 
of spatially separated packets, and of the energy production rate for any number of free particles: 
(N ■ these reduce to the CSL rates if the particle's Compton wavelength is small compared to the model's 

. distance parameter. One motivation for QRCSL is the realization that previous relativistic models 

' entail excitation of nuclear states which exceeds that of experiment, whereas QRCSL does not: an 

Lj^ | example is given involving quadrupole excitation of the 74 Ge nucleus. 
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CNj \ The Continuous Spontaneous Localization (CSL) model0, is, currently, the only modification of Schrodinger's 
equation which satisfactorily describes both standard quantum physics and the macroscopic world we see around 
us. When a superposition of macroscopic spatially distinguishable states appears, the wave function dynamically 
. undergoes rapid collapse toward one such state. 

Since the world is locally Lorentz invariant, it is natural to try to make a relativistic collapse model. I have been 
working at this for over a decade and, although the models constructed have certain interesting features, each has 
■ certain interesting flaws. 

The major problem with the first models^, 0, is that, even to lowest order in the collapse rate parameter A, 
there is an infinite energy production rate per volume from the vacuum. The reason is as follows. In nonrelativistic 

j— i CSL, a randomly fluctuating classical field w(x,t) interacts with the mass-density operator (coupling constant A). 

• • , Collapse narrows wavepackets, resulting in a small rate of energy increase of particles (the energy is supplied by 
w(x, t) 6]) which is, at present, below experimental observation^. However, in these relativistic models, each vacuum 
energy-momentum mode, which may be regarded as the ground state of a harmonic oscillator, is likewise narrowed by 
the collapse interaction. Thus each vacuum mode is no longer just the ground state but is a superposition of ground 

?3 ! an d excited states, i.e., there is a small rate of creation of particles with that energy-momentum. Since there are an 
infinite number of modes, than infinite energy production rate per volume of particles from the vacuum. 

In nonrelativistic CSL and in these first relativistic models, w(x,t) is white gaussian noise, i.e., it contains all 
wavelengths and frequencies in equal amounts. In standard quantum theory, interaction of an operator with classical 
noise of frequency u) results in excitation of the quantum system with energy change huj. That is also the case for 
collapse models. 

I therefore considered models where the noise is non- white gaussian 8, 9] (which entails a non-Markovian state vector 
evolution) and showed that, indeed, suppression of the frequency u> in the inverse of the spectrum of w(x, t) suppresses 
energy hu> excitation. This led to a relativisitic model with a tachyonic (momentum-energy relation k 2 - fc° 2 = -fi 2 ) 
inverse noise spectrum|lOj, which does not excite the vacuum to lowest order in A since there is no mode of the vacuum 
which has tachyonic energy-momentum. 

However, to order A 2 , once again the ugly spectre(um) of infinite energy production from the vacuum appears. The 
culprit is that, in this order, there is an internal particle line in the Feynman diagram describing vacuum excitation. 
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The associated particle propagator, (k 2 — fc° 2 + M 2 — ie)" 1 , like the white noise spectrum, possesses all wavelengths 
and energies, which it uses to convert the tachyonic energy-momentum to vacuum production of a pair of particles. 

If the particle propagator were on-shell, <5(k 2 — k° 2 + M 2 ), then the vacuum excitation disappears to this and 
all orders. This may be achieved by removing the time-ordering operation from the state-vector evolution operator, 
resulting in a finite relativistically invariant collapse model (RCSL), but at a cost. First, the model is expressed 
in the interaction picture: with time-ordering one can reconstruct the Schrodingcr picture (multiply the interaction 
picture state vector by exp —iHt, where H is the Hamiltonian) but, without time-ordering, this connection is severed. 
Second, the model is nonlocal in an unusual sense since, with time ordering, each Feynman diagram describes a 
series of forward- in-time evolutions but, without time-ordering, these alternate with backward- in-time evolutions. 
While the time-ordering evolution keeps an evolving particle rigorously within the light cone of its initial spatial 
state, the non-time-ordering evolution does allow the particle to go out of the light cone in high enough order, albeit 
with a small probability. Third, and most decisively, as shown in Section V and Appendix A, to lowest order in A 
(where time-ordering plays no role), the tachyonic-based theory predicts too great an excitation of nuclear states: the 
"spontaneous" quadrupole excitation from the ground state (0 + ) of 74 Ge to the first excited state (2 + ) at .596MeV 
greatly exceeds the experimentally observed rate. To my knowledge, this is the first situation where experiment has 
ruled out a collapse model. 

This paper takes another tack. I have been unable to to construct a collapse model simultaneously satisfying 1) 
relativistic invariance, 2) locality, 3) Markovian evolution, 4) no vacuum excitation, 5) consistency with experiment. 
Nonrelativistic CSL satisfies all but 1), replacing it with galilean invariance. Relativistic collapse models in refs0,0>@ 
do not satisfy 4), 5) and RCSL does not satisfy 2), 3) and 5). (See references [Tit IT^ ] for a similar assessment of recent 
interesting models by Rimini and Nicrosini and by Tumulka.) This paper contains a model which satisfies 3), 4), 5), 
which satisfies 1) and 2) approximately, and which reduces to CSL in the nonrelativistic limit. 

II. NONRELATIVISTIC CSL 

Nonrelativistic CSL is characterized by two equations, the state vector evolution equation and the probability rule. 
The evolution equation in the "collapse interaction picture" (where the operators are Heisenberg operators, and 
the state vector only changes with time due to collapse) is 

|V, t) w = Te-^ 1 /o ^xMx, t )-2AA(x,t)] 2 Q) (1) 

(T is the time-ordering operator). In Eq.(l), A(x, t) is an operator essentially proportional to the mass of particles 
in a spherical volume of radius a, and can be written in various ways: 

A(x,0) ee Y,^™ 2 ) m e 2 - la2v2 Zl(x)U*), (2a) 
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dbe-( 2a2 )" lb2 4(x + b)£ a (x + b) (2b) 



= E W (iJ) I ^P'^- 2 " a2(P '- p)2 e- 4(p '- p) x 4(P> Q (p)- (2c) 

where A(x, t) — exp(zi/t)A(x, 0) exp(— iHt), with H the complete Hamiltonian for the interacting particles. In 
Eqs.(2), 

£t (x )^(2^r3/ 2 ! dpe-^aUp) (3) 

is the creation operator for a particle of type a (electron, proton, neutron,...) at the position x (so £J,(x)£q,(x) is the 
number density operator for particles of type a), a„(p) is its momentum p creation operator, M a is its mass and M is 
the proton mass. The values of the parameters A and a which characterize the model are generally chosen to be those 
given in the seminal collapse model of Ghirardi, Rimini and Weber (GRW) 01, A = 10 _16 sec _1 and a = 10~ 5 cm, 
but there is a range of parameter values allowed by experiment 0] . 

The probability rule gives the probability density that the classical noise field w(x,t) occurs in nature: 

V(w,t) EE w (lP,t\lp,t) w . (4) 

That is, Eq.(l) does not preserve statevector norm, and Eq.(4) says that state vectors of large norm are most likely to 
occur. It follows from Eqs.(l), (4) that J DwV(w,t) = 1, where Dw = Y[ x t dw{~x.,t){2ir\/ 'dndt)^ 1 / 2 is the functional 
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integration volume element (in doing the integrals, x, t are discretized, with w(x, t) regarded as an independent 
variable for each (x, t). 

It is readily shown that Eqs.(l), (4) entail that a state vector, describing a macroscopic object in a superposition of 
places, rapidly evolves toward one of the states in the superposition with probability equal to the squared magnitude 
of its coefficient in the superposition. Essentially, a state which survives when all others have collapsed is one for 
which the time average of w(x,t) equals the state's time average of 2XA(x,t) at each x (all other behaviors have 
vanishing probability) . 

For calculations of physical effects, it is easiest to utilize the density matrix which describes the ensemble of 
evolutions: 

p(t)= Dw————V(w,t) 5a 

j w{i>,t\ip,t)w 

= Te~ {X/2) - f o dtdx{A L (x,t)-A R (x,t)] 2 p ^Q^ ^ 

Eq.(5b) follows from putting Eqs.(l), (4) into (5a). In Eq.(5b), when the exponential is expanded in a power series, 
Al (Ar) operates to the left (right) of p(0), and T time-orders (time-reverse orders) the operators at the left (right). 



III. QUASIRELATIVISTIC MODEL FOR NONINTERACTING PARTICLES 

In what follows, for simplicity, only one type of particle, a "bosonic nucleon" of mass M shall be considered: 
the results are trivially extendable to fermions and bosons of many types. The relativistic generalization of the 
nonrelativistic creation operator £'(x, t) for noninteracting particles (Eq.(3) where expiEt is stuck into the integral, 
with E — p 2 /2M) is the negative frequency field operator 

<T(x,t) = (2ir)-V 2 f dv{M/Ef^e- ip - x (J{v) (6) 



[p ■ x = p • x — Et, where E = (p 2 + M 2 ) 1 / 2 ) which, like the local field <p{x) = <fi (x) + <j) + (x), transforms like a scalar 
under Lorentz transformations. In Eq.(6) we choose the commutation relation [a(p), a^(p')] = S(p — p') so a(p)>/E 
transforms like a scalar. It is clear that </) _ (x, t) reduces to £'(x, t) in the nonrelativistic (c — > oo) limit (except for 
an additional factor expiMt, which cancels out in <fi~(j) + ). 

In this quasirelativisitic CSL (QRCSL) model, the evolution equations are also (1) for the state vector, (5) for the 
density matrix, and (3) for the probability rule. But, for QRCSL, A(x, t) is defined using an approach from nonlocal 
relativistic quantum field theory [l^, and can be written in various ways parallel to Eqs.(2): 

A(x) = {ATra 2 f/ A e 2 ^ a2n (t>-{x)(t> + {x) (7a) 

dbdb e~ {2a ' 2rl[b2+b « ] (j)- (x + b, t + i6 o )0+(x + b, t + ib ) (7b) 



2 1 /2( 7ra 2)5/4 

= (iJ) 3/ Y ^P-^e- 2 -V-P) 2 aV^'-ri-V(p>(p) (7c) 

( a ^ 2 -d 2 ). 

It is apparent from Eqs.(7a,c) that A(x) is a Lorentz scalar. It is also worth noting that the exponent (p 1 — p) 2 in 
(7c) is spacelike (i.e., positive: in the reference frame where p' = 0, (p 1 — p) 2 — 2M(E — M)). It is easy to see from 
Eqs.(7) that, in the c — > oo limit, A(x) reduces to the nonrelativistic A(x,t) (Eqs.(2) with time behavior added). 

The equal time commutator, utilizing Eq.(7b), is 

[ J 4(x,i), J 4(x',t)] - (2^/V)- 1 J dhdbadh'db^e-^^+^e-^ 2 ^ 1 ^' 2 ^ 

[(f>-{x + b,t + ib )(f> + (x' +b',t + ib' ) - (j)-{x +b',t + ib' )(f> + (x + b,t + ib )} 

% 

with 



+ (x + b,^ + ^6o),^(x' +b',t + ib' )] (8) 



[4> + (xU-{x')] = -^ J ^ c *-C*-» '). ( 9 ) 



4 



Eq. (9), with the arguments appropriate to (8), is 

, +( . +b , t+it0) , r( , + ,, t+a ^ ^^^^g . (10) 

where K\ is the Bessel Function. Now, b,b',6 ,^o m Eq.(8) have gaussian distributions with spread a. For most 
of their volume of integration where the gaussians are large, the argument of K\ in (10) is quite large, of order 
Ma = 10~ 5 cm/10~ 14 cm= 10 9 . Since for large argument, K\(z) — > (■n/2z) 1 / 2 exp — z, in Eq.(8) where the gaussians 
are large, the factor K\ is, for the most part, quite small. It is in this sense that the commutator (8) "almost" 
vanishes, making the time-ordering operation "almost" frame- independent and the model "quasi" -relativistic. 
In calculations of the density matrix using Feynman diagrams, internal particle lines are represented by 

(0\T[<t,+ (x),<t>-&)]\0) = (*-*')^3 / ^e*-<*-*'> (11a) 

_ [ di P 1 ,»•(»-»') fl lb) 

" (27r) 4 i E p*-E + ie { ' 

2Mi [ A H[(/-y] e ,M'), (llc) 



(2tt) 4 J V 2 -p 2 -m 2 +ie 

Eq.(lla) follows from Eq.(9): in (11a), p = (p,E) while in (llb,c), p = (p,f>°). 

Eqs.(ll) show the lack of Lorentz invariance. However, with (x— x', t—t') replaced by (x+b— x' — b', t+ibo — t' — ib' ) 
in (11a), the integrands of (11a, b,c) acquire the factor expi[p-(h—h')+E(bo — b' )]: for spacelike (x—x'), the propagator 
equals Q(t— t') multiplied by Eq.(10) with argument [M(x + b x' b') 2 — (t+ ib a —t' — ifr^) 2 ] 1 / 2 ), and the previous 
discussion of how this "almost" vanishes applies. 

One might consider writing the propagator in each Feynman diagram, and thus the whole density matrix, as the 
sum of a relativistic piece plus a non-relativistic correction. There are various ways to achieve this. One might split 
the propagator in space-time into its expression within the forward light-cone and zero elsewhere, with the correction 
as the spacelike remainder. Eq.(llc) shows another split, in momentum space, with the relativisitic part equal to 
the usual propagator (recall that our definition (6) of <f>~ is (2M) 1 / 2 times the usual definition) plus a part which is 
relatively small for p a w E. 



IV. CALCULATIONS 



The density matrix evolution equation and its perturbation series follow from Eq.(5b): 
dp(t) 



= -(A/2) J dx[A(x,t),[i4(x,t),p(t)]] (12a) 



dt 

Pit) = V^l^ f dx n ... ( dx 1 T[A(x n ),...[A(x 1 ),p(t)}...}. (12b) 

In this section, Eq.(12a) is used to calculate the rate for the wave function of a single free particle, in a superposition 
of two widely separated packets, to collapse to one of the packets. Next, the rate of energy increase for N free particles 
is found. The following section discusses the formalism of QRCSL when particles are interacting. Finally, this result 
is used to calculate, to first order in A, the quadrupole excitation rate of the 74 Ge nucleus from its ground state to its 
first excited state. The result is compared with the present experimental upper limit on the rate of "spontaneous" 
excitation in Ge. 



A. Collapse Rate for a Single Free Particle 



Consider a single particle initially in a superposition \ip, 0) = a\L) + (3\R), where \L) and \R) are widely separated 
wavepackets, so far apart that, to high accuracy, their regions of support do not overlap over the time interval t. Let 
|xl) (|x#)) be a position eigenstate within the left (right) region of support. The off-diagonal element of the density 



■5 



matrix, using Eqs.(12a) and (7c), satisfies: 

<HW M i dx I 



d(x ^' f )|Xfl) = -(A/2)f4V /2 ^ 2 / f ^1^^^ 

dt y 1 J \Air 3 J J J y/El y/Ei y/Ei y/Ei 

2~ 1 a 2 [(p 1 -(p 2 ) 2 + (p3-(p4.) 2 ] p -i(Pl-P2+P3-pi)-x 



{ 5 (P2-P3)(^)3 / dz{e^- 1 ^(z\p(t)\x R ) +e^-^«(x L \p(t)\z}} 

~ 2 J^Y I dz J rfz ' e ' pl ' XI '^ P2 ' z+iP3 ' z ' _iP4 ' Xi? ( z l/ 9 Wl z ')|- 



(13) 



First, it can be seen that the last term in the curly brackets of Eq.(13) (arising from the term 

— A J dxA(x,t)p(t)A(x,t) in Eq. (12a)) may be neglected. Consider the integral over p 1; p 2 which appears in 
this term: 

/(X£,X,Z;t) = / ^EL-^ =e - 2_la2 (Pl-P2) 2 e - 4 (pi-p2+P3-p4)-a ;e ip 1 -X L - t p 2 .Z ^ 14ft J 

» ^ e -( 2a2 ) _1 ( x - x ^ 2 ,5(z-x L ), (14b) 

where w in (14b) means that we have set t = and E\ = E 2 = M. Even without the approximation, for any t, if any 
one argument (x or z) of / lies in L and another in R, then / w 0. Now, the term under consideration has the form 

~ [ dxdzdz' f(x L ,x,z;t)f*(x R ,x,z';t){z\p(t)\z') (15a) 

= J dxdzdz' f(x L ,x,z L ;t)f*(x R ,x,z' R ;t)(z L \p(t)\z' R ) ^0. (15b) 

In (15b), z (z') is restricted to z^ (z' R ) since, otherwise, / (/*) would vanish, and the result vanishes because either 
/ or /* vanishes for every x . 

In the remaining terms of (13), first perform the integral over x, followed by the integrals over p4, which results 

in 



d(x L \p(t)\x R )_ = _( A / 2 ) 



dt 



J E\ Ei 

J dz[e l P 1 -( Xi - z )(z| (0 (t)|x K )+e- 4 P 1 -( x "- z )(x L |^)|z)]. (16) 
Next, the integral over p 2 can be performed: 

f ^l e -a 2 ( Pl -p 2 r = ( 27r /a 2 )e 2 - aM V 1 (2a 2 M 2 ) ) (17) 

J &2 

followed by the integral over pi: 

/ ^e* 1 "^-) = 4ttM^ 1 (M|x l - z|)/|x L - z|. (18) 
The evolution equation now reads: 

' i(Xtl l" |X ° > = -{XaM'ftr'IV'^K^aHf-) 

f ^iMW lj) (, Wl) |,„) + gh ■-■'' { W |p(0|.)]. (19) 

J | X L — Z| |X^J — Z| 



Eq. (19) is exact. 

Now, specialize to the case where a is much larger than the particle Compton wavelength M^ 1 so 

e 2a 2 M 2 K ^ 2a 2 M 2 ) M ( 2a M)-V 1 / 2 . (20) 
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Further, if the wave packets change slowly on the distance scale of M -1 then, in the integrand of Eq.(19), K\{M\x — 
z|)/|x — z| may be well approximated by a delta-function, whose numerical coefficient may be found by integrating 
Eq.(18) over z: 

ifi(M|x - z|)/|x - z| w (2ir 2 /M 2 )6{x - z) (21) 
The result of inserting (20), (21) into (19) gives the result: 

d{XL m\ XR ) = _ x{XLlp{t)lXRh (22) 

i.e., rate of decay of the off-diagonal density matrix element is A, the same rate as for nonrelativistic CSL. 



B. Energy Creation Rate For Free Particles 



Calculation of the energy creation rate for n free particles begins by multiplying Eq.(12a) by the free Hamiltonian 
H and taking the trace: 



dH_ 
~dt 



A 
2 

\i 

y 



J dxTr{[A(x),[A(x),H ]}p(t)} 
dxTr{[A(x),A(x)]p(t)}, 



(23a) 
(23b) 



where H(t) = Tr{Hop(t)}. Using the expression (7c) for A(x), the commutator can be evaluated. Then, the integral 
over x can be performed and, using the resulting delta-functions of momentum differences, two of the momentum 
integrals can be performed, resulting in 



dH(t) _ XM 2 a 3 f dpi dp 2 



dt 
AM 2 " 3 



J E 1 



'-(E 2 -E 1 )e 



-a (pi-p 2 ) 



2tt 3 / 2 J Ei E 2 

Tr{[at( Pl )a(pi) - a\p 2 )a(p 2 )}p{t)} 

^Tr{[/|l a t (pi)a(pi) J dp2e -^- P >f- N J ^e-^-^W)} 

In Eq.(24b), N = J (i P ia' l '(pi)a( P i) is the number-of-particles operator. 

The last integral in Eq.(24b) over p 2 is given in Eq.(17). The first integral over p 2 is 

J dp 2 e- a2 ^- p ^ 2 = (2ir/a 2 )E ie 2a2M2 [K Q (2a 2 M 2 ) + (a 2 M 2 )" 1 K^a 2 M 2 )] 



(24a) 
(24b) 



(25) 



Because (25) is proportional to E\, the first integral in Eq.(24b) over P i is ~ N. Since Tr{Np(t)} = nTrp(t) = n, 
one obtains the exact result that H(t) increases linearly with time: 



= AnM 2 a 27r- 1 / 2 e 2a2M2 [^ (2a 2 M 2 )-X 1 (2a 2 M 2 )(l-(a 2 M 2 )- 1 ) 



In the case where aM » 1, using the series expansions 



Eq.(26) becomes: 



dH(t) 
dt 



3An 



1 + l 



4Ma 2 ' 



(26) 



(27) 



(28) 



the same energy creation rate as for nonrelativistic CSL. 
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The derivative of Eq.(l), 



V. QRCSL FOR INTERACTING PARTICLES 



d\i>,t) w /dt = -(AX)- 1 J dx[w(x) - 2\A(x)] 2 \tP,t) w (29) 



where A(x) is given by Eqs.(7), is the QRCSL evolution equation for noninteracting particles in the "collapse interac- 
tion picture," where the field operators evolve freely and the state vector changes only due to the collapse evolution. 
As usual, the Schrdingcr picture statevector is = ex P — iH t\ip, t) w , and its evolution equation follows from 

Eq.(29): 

d\ip,ty w /dt = -tH \^,ty w -(4X)- 1 J d x [w(x)-2XA( x ,o)} 2 \^ 7 ty w , (30) 

where A(x, 0) is given by Eqs.(7) with t = (so, in particular, the operators in (7b) still have ibo as time arguments). 
In the Schrodinger picture, the operators do not evolve, and the state vector changes with time due to the free 
evolution as well as due to the collapse evolution. 

As usual, to add interaction, one replaces H by H = H + V in Eq.(30). Transforming back to the interaction 
picture gives: 

d\^,t) w /dt = -iV(t)\i>,t) w - (4A)- 1 j dx[w(x) - 2XA(x)] 2 \ij,t) w , (31) 

where V (t) = exp(iHot)V exp — (iH t) is a Lorentz scalar, the four-integral of a local scalar density. However, although 
Eq.(31) is form covariant, it is not Lorentz invariant because, not only doesn't A commute with itself at space-like 
separation, it usually will not commute with the local scalar density at space-like separation either. However, such 
a local density is constructed from 4> = 4> + + <t>~ > so its commutator with A falls off exponentially with space- like 
separation as in (10). In this sense this interacting QRCSL model is also "quasi-relativistic." In this interaction 
picture, the quantum fields evolve freely and the state vector evolves due to the interaction and the collapse. 

For some calculations, it is useful to work in the collapse-interaction picture, where the fields arc Heiscnbcrg fields, 
evolving according to the interacting quantum field theory and the state-vector only changes with time due to the 
collapse evolution. Using (31) to go to the Schrodinger picture and, as usual, defining the collapse-interaction picture 
statevector as \ip,t)' w — exp(iHt)\ip 7 1)^, one obtains 

d\^t)'Jdt = -(4A)" 1 / dx[w(x) -2\A'(x)] 2 \^ 1 t)' w , (32) 

where A' (x) = cxp(iHt)A(x, 0) cxp —(iHt). It is worth emphasizing that A' (x) is not given by Eq.(7b) with operators 
^(x + b, t + ibo) replaced by exp[iH(t + ifeo)]0 ± (x + b, 0) exp — [iH(t + ibo)] but, rather, is composed of operators 
exp(i_fft)0 ± (x + b, ibo) exp— (iHt), according to this prescription. 



A. Bound State Excitation to Lowest Order 



The density matrix evolution equation which follows from Eq.(32) is Eq.(5b) with A replaced by A'. To lowest 
order in A, this is 

P (t) = p(0) - (A/2) / d^xL4' 2 (x, t) P (0) + P (0)A' 2 (x, t) - 2A'(x, t))p(0)A'(x, t)}. (33) 
Jo 

Time ordering plays no role to lowest order, so the expression (33) is Lorentz invariant. Take p(0) = \Ei)(Ei\, where 
\Ei) is a bound A^-particle state that is an energy eigenstate with eigenvalue Ei, and is also an eigenstate with 
eigenvalue of the center of mass operator Q = J2n=i X n /A (X„ is the position operator of the nth particle). It is 
desired to calculate the probability that the system is excited to the bound energy eigenstate \Ef) (H\Ef) = Ej\Ef) 
and Q\Ef) = 0). Because (E^Ef) = 0, only the last term in (33) contributes: 



(E f \p(t)\E f ) = A f dtdx{E f \A'(x 1 t)\E t ){E t \A'(x 1 t)\E f ). 
Jo 



(34) 
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Moreover, since (E t\A' (x, t)\Ei) = expi(Ef — Ei) (E f\A(x, 0)\Ei) , the integrand in (34) is time independent, so the 
excitation rate T = d(E f\p{t)\E f) / dt is constant: 



r = \Jdx\(E f \A(*,o)\E l )\ 2 . (35) 

At this point, assume that the particles in the initial and final states move nonrelativistically, (p/mc) 2 << 1, so 
that, in the expression (7c) for A(x), one can make the approximations E m E' s» M and (p — p 1 ) 2 ~ (p — p') 2 - Then 
A(x, 0) in Eq.(35) becomes the non-relativistic expression (2). For completeness, the analysis leading to Eq.(37) shall 
be given here, instead of just auoted [l4L Il6[ . Start by noting that J <fx'.F(x')£t(x')£(x') = J2n=i F(X n ) where F is 
an arbitrary function. Use of the form (2b) for A(x, 0) in Eq.(35) results in: 



= A(™ 2 )- 3/2 J dx\(E f \ J dx' e ~( 2a > < x - x ') £t( x ')£( x ')|£.)|2 (36a) 
r N 

= A(™ 2 )' 3 / 2 / dxK^l^e-^" ) ( x - x «)> i )| 2 (36b) 

^ 71=1 

N 

= A J2 {E f \{e-^- 1 ^-^ 2 \E,){E l \}\E f ) (36c) 



n,m— 1 



Next, expand the expression in (36c) in a power series in (size of bound state/a) 2 . Then, because of the orthogonality 
of the initial and final states and because Q|-Ej,/) = 0, the first nonvanishing term of (36c) is of order a~ 4 : 



T = A(2a)~ 



N N 3 



\(E f \^2x 2 n \E i )\ 2 + 2 £ ^(^ix^x^i^^ix^i^i 

n— 1 n— l,m— 1 — l 



(37) 



B. Excitation of Ge Nucleus 



Now, Eq.(37) is to be applied to collapse- induced spontaneous excitation of a proton from the ground state (0 + ) 
of a 74 Ge nucleus (the largest percentage isotope-36.5%-in naturally occurring Ge) to its first excited state (2+) 
at .596MeV^3- The experimental upper limit on spontaneous emission of .596MeV gammas in Ge, obtained by 
observing the radiation from an isolated slab of Ge for a long time, is w .03 counts/kg-day (in a 2 MeV bin) The 
analysis for spinless particles of one mass given above can be applied to this case because the proton and neutron 
have almost the same mass and there is no spin-flip involved in this transition. 

One readily finds from (37) the expression for the quadrupole excitation rate: 

2 z 

T = (7r/15)(A/a 4 ) £ |<2+, m'\ ]T X 2 r 2m (6„, $„)|0+)| 2 , (38) 

m,m' — — 2 n—1 

where ra , <fr ra are angle operators for the nth particle and Yim is a spherical harmonic. 
Now, the lifetime r of the 2 + state is given by the expression [l9j 

2 z 

r- 1 = (47r/3-5 3 )cfc 5 (e 2 Ac) £ |(2+, m'| £ X 2 F 2m (e„, $„)|0+)| 2 , (39) 

m,m' — — 2 n—1 

where k ps 3.2 • 10 10 cm _1 is the .596 MeV photon wavenumber. From Eqs. (38), (39) is obtained 

Fqrcsl = (5/2) 2 X[(ak) i (e 2 /hc)kcT}- 1 w 5 ■ 10" 16 counts/kg-day, (40) 

with use of the numbers r = 17.9psec as the experimental lifetime of the state, w 8.3 • 10 24 as the number of nucleii/kg 
of Ge (so there are w 3.0 • 10 24 74 Ge nucleii/kg), and 8.6 • 10 4 sec/day. Tqrcsl is well below the experimental upper 
limit of 3 • 10~ 2 counts/kg-day. 

This contrasts with the situation for the relativistic collapse model RCSL. From Eqs.(A4) and (39) is obtained 

Trcsl = (5/3 7 r 2 )Aa[(e 2 /?ic)cT]~ 1 « 5 • 10 10 counts/kg-day, (41) 
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which far exceeds the experimental upper limit. These calculations were performed assuming the GRW values for A, 
a but, for most of the range of these parameters allowed by other considerations ^4|, Trcsl is still excessive. The 
reason for the difference in excitation between QRCSL and RCSL is that, in the former, just as in nonrelativistic CSL, 
collapse narrows the excited particle's wavefunction to a whereas, in the latter, collapse narrows the wavefunction to 

VI. CONCLUDING REMARKS 

Because all previous CSL-type relativistic collapse models except RCSL are untenable since they produce infinite 
energy /sec- vol from the vacuum, and RCSL produces too much nuclear excitation, the QRCSL model has been 
suggested. It has form-invariant equations, but it fails to be relativistic because its Lorentz invariant operators 
do not commute at space-time separation. However, since these operators "almost" commute, I believe that such 
quasi-relativistic behavior is worth consideration, as a close and experimentally testable variant of special relativity 
combined with a description of collapse that is as close as could be expected to nonrelativistic CSL. 

However, the model is, after all, described in a preferred frame, the one where the time-ordering operation is defined. 
One might tentatively identify the preferred frame with the local co-moving frame of the universe[2(|. Exploration of 
the extent of violation of Lorentz invariance for various hypothetical situations is certainly of interest. Since QRCSL's 
slow speed limit is CSL which, so far, has defied experimental refutation, it may be worthwhile to examine schemes 
whereby detectors move at high speeds. One may also examine whether the frame dependent, although non-detectable, 
wave packet collapse locales and times (e.g., in EPR-type situations) in relativistic collapse models are similar to those 
of QRCSL, or if the preferred frame's wave packet collapse locales and times might, in some sense, be preferred 

APPENDIX A: BOUND STATE EXCITATION IN RCSL 

In RCSL, the only finite relativistic collapse model extant, the expression for the excitation probability comparable 
to Eq.(34) (but here taken in the nonrelativistic limit) is 

(E f \p(t)\E f ) = 4Aa f dxdx'G{x- x'){E f \(}{x)^x)\E l ){E l \e{x l )^x l )\Ef). (Al) 
Jo 

In Eq.(Al), G(x — x') = (27r)~ 4 J dpexpip ■ (x — x')8{p 2 — a -2 ), i.e., this is a non-Markovian model with a tachyonic 
noise spectrum whose "tachyon mass" is a" 1 w 2eV. The operator £'(x) = expiifi£t(x, 0) exp — iHt, with £t(x, 0) 
given by Eq.(3), is the Heisenberg creation operator. Use this first, to pull out the time dependence from the matrix 
elements in (Al), letting H act on the energy eigenstates, and then perform the time integrals with the result 

(E f \p(t)\E f ) = 4(27r)~ 4 Aa f dxdx'dpe lp{x ~ x,) S(p 2 - a" 2 ) 

Jo 

•{ sm[(E f - Ei-p )t/2]/[(E f - E t -p°)/2] j {E f \e^)^m){E^^)^')\E f ). (A2) 

For large t, [sin(crf)/a] 2 « tird(a). Then, using this delta function to perform the integral over p , and utilizing 
J dx'F(x')£t(x')£(x') = 5Z n =l F(yZn) as was done in obtaining Eq.(36), we obtain 

N 

T = Ao2- 1 tt- 3 / dp5[p 2 - (E f - E % ) 2 - a' 2 ] £ {E f \{e^^-^\E l ){E t \}\E } ) (A3a) 

n,m— 1 

N 

= \air- 2 (Ef\{[MQ^-^R\)/\^L-XmR\]\BH){Ei\}\Ef). (A3b) 

n,m— 1 

In Eq.(A3b), k — yj [Ef — Ei) 2 + a~ 2 w Ef — Ei (if Ef — Ei >> a -1 ) is the wavenumber of a photon making the 
transition from the excited state to the ground state. 

Eq.(A3b) may be compared to Eq. (36c). The gaussian with width a there, is replaced by the sinfcz/z form with 
width fc" 1 here. An expansion in powers of a -1 there is replaced by an expansion in powers of k here. It is because 



10 

k >> a^ 1 that this RCSL model produces a much larger excitation rate than QRCSL. The first non-vanishing term 
in Eq.(A3b) is 



N N 



\{E f \j2K\E<)\ 2 + 2 Y.^ E f\^m)m^ i m ^ m \E f \ 

n—l n— l,m— 1 i,j=l 



r = 2Aa/c 5 (5!7r 2 )- 1 

This is identical to Eq.(37) except for the numerical factor, and is used in Eq.(41) 



(A4) 
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